Solution of the two identical ion Penning trap final state 
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We have derived a closed form analytic expression for the asymptotic motion of a pair of identical 
ions in a high precision Penning trap. The analytic solution includes the effects of special relativity 
and the Coulomb interaction between the ions. The existence and physical relevance of such a final 
state is supported by a confluence of theoretical, experimental and numerical evidence. 



PACS numbers; 32.80.Pj, 02.20+b, 33.80.Ps 

High precision Penning traps are ideal for studying phys- 
ical characteristics of individual ions. These traps, as 
described for example in Ref.[l], have magnetic fields 
that over the trajectories of the ions vary by less than a 
part per billion. In consequence, the motional frequency 
linewidths can be made so narrow that effects of special 
relativity are readily apparent even at these relatively low 
velocities^. 

To remove systematic effects it is often desirable to fill 
the trap with two ions and much is known about the re- 
sulting frequency perturbations caused by the Coulomb 
interaction between dissimilar ions^. The situation with 
two identical ions has also been extensively studied much 
(see Ref.[4,5] and references therein). The solution and 
approach that we describe here are rather different than 
those references however, since they include the electric 
trap field but ignore relativistic mass increase. Includ- 
ing this effect of special relativity may be crucial for un- 
derstanding the observation^ of cyclotron mode-locking 
between identical ions (see also Ref.[7]). 

We present details of an analytical model of two iden- 
tical ions in a high precision Penning trap. The model 
is asymptotically solvable in terms of elliptic functions. 
This solution is, in practical terms for protons and heav- 
ier ions, a generic final state of two identical ions in a 
precision Penning trap. 

We begin with a symmetry argument detailing what 
is special about the two identical ion system and then 
we introduce and solve the model. For two dissimilar 
ions the center of charge is different than the center of 
mass. The motion of the center of charge causes currents 
to run in the detection circuit and in the walls of the 
trap itself causing a force to act back on the ions. This 
retarding force acts on the center of charge and so if the 
center of charge is different than the center of mass these 



damping forces act always on a mixture of the center of 
mass motion and the relative motions of the ion pair. 

This is not the case for identical ions in the trap. In 
that case the center of mass and the center of charge 
are the same and so the retarding force acts only on the 
center of mass motion. Thus, the relative motion of the 
ions is relatively undamped, being subject only to the 
weaker quadrupolar damping (which is associated with 
timescales generally longer than typical experiments). In 
this sense we speak of this final state of the two identical 
ion system as a decoupled, or, dark state. 

One way to understand the existence of this cyclotron 
dark state is with a symmetry argument. Neglect dissi- 
pation, relativity and interaction and consider the Pois- 
son algebra of two ions moving in a horizontal plane (we 
shall describe why this is relevant to experiment later) 
in a uniform perpendicular magnetic field. The Hamilto- 
nian is proportional to Ti. — pi + P2 + ct{P3 + p1), where 
a = ma/mi, is the mass ratio and pi^2 (resp. ^3,4) are the 
canonical momenta of particle a (resp. particle b). For 
a ^ 1 the subalgebra commuting with Ti. is so(2) x so(2) 
whereas if a = 1 the algebra is so(2) x so(3). The fact 
that there are additional commuting generators in the 
equal mass case indicates that there is a flat direction in 
the dynamics of that case, corresponding to degeneracy 
between cyclotron dark states of different total angular 
momentum. 

There is a straightforward geometrical way of under- 
standing the special qualities of the two identical ion 
Penning trap. Again consider the ions confined to a 
plane perpendicular to the magnetic field and ignore tem- 
porarily the effects of relativity and interaction. The to- 
tal angular momentum of the two ion system is L = 



Pi 



P3 + Pi (note independent of the mass ratio 



a). Now, turning on relativity and interactions pertur- 



batively, we learn that the motion is essentially restricted 
to the intersection of iso- Ti and iso- L surfaces. A generic 
intersection of these surfaces in R'' for the a 7^ 1 case is 
a two-dimensional torus (and so has an isometry group 
so(2) X so(2)) whereas when a = 1 the intersection is 
not generic, but is the whole . Although the isome- 
try group of , being so(4), is isomorphic to so(3) x 
so(3) the physically relevant isometry group is that which 
preserves not only the geometry but also the underly- 
ing Poisson structure, which is sp{A) in this case. The 
canonical intersection^ in the group of matrices GL(4) of 
so(4) and sp(4) is the algebra m(2), which is isomorphic 
to so(2) X so(3), which again is the enhanced symmetry 
discussed above. We note that both the geometrical and 
algebraic picture can be easily generalized to the case of 
A*" identical ions^. 

Having described the symmetry properties unique to 
two identical ions in a Penning trap, we now introduce 
the interacting model by starting with the following three 
assumptions. 

1) The ions are very near the center of the trap, and 
ignore effects due to the spatial gradient of the electro- 
static fields of the trap (that is, we completely ignore the 
trap magnetron motion). The cyclotron frequency shifts 
in an isolated ion's cyclotron motion is entirely due to 
relativistic effects. 

2) the ions are mode locked already in the trap's axial 
drive and so their motions may be thought of as being 
confined to a plane®'^. 

3) The energy loss mechanism is entirely due to the 
dissipation of image charge currents induced in the 
trap/detection system, and thus couple only to the center 
of mass of the ion pair) . 

Under these assumptions, the equations of motion for 
the ion pair are the formidable looking non-linear coupled 
differential equations; 



fi + ujq{\ - fi)z X f 1 + 7A„„ - 



r2 + c^o(l - f2)z X r 2 + 7A„„ + 



= 



= 



(1) 



(2) 



where Xcm — (n. + r2)/2, and i? = ri — r 2 and where 

I ^ 1 2 

/i — '-^jr is just the ratio of the kinetic energy to the 
rest mass-energy of ion 1 (similar expression for /2 is in 
terms of the kinetic energy of the second particle). This 
term, due entirely to special relativistic mass increase, 
causes the cyclotron frequency to depend on the kinetic 
energy of the ion(s). 

We add and subtract Eq. (|l]) and Eq. (||) to rewrite 
them in terms of the center of mass co-ordinate Xcm and 
the relative coordinate R, 



Xcm + "^^Xcm + ^^o(l 



/l + ./2 



)z X Xc 



R + uJo{l 



(4) 



Let V = Xcm be a symbol for the center of mass velocity. 
As expected, only the center of mass velocity enters into 
the equations. Confined as they are to the same vertical 
plane, this becomes a six-dimensional (phase-space) sys- 
tem. Let U — R. In these variables, the combinations 

/i-/2-^andA + /2 = ^^. 

As per earlier discussion, from Eq. (||) and Eq. (^, it 
is clear that the center of mass motion is damped but 
the relative motion is not. Thus, after sufficient time, it 
is consistent to assume that the center of mass motion 
damps out completely, that is, V —^ 0. The coupling 
term between the R motion and the V (center of mass) 
motion is through the term proportional to /i — /2 (itself 
proportional to V), and so Eq. (^ and Eq. (^ quickly 
decouple as V 0. 

The resulting motion can be treated perturbatively in 
small V. To find the zcroth order term we ignore the cou- 
pling term completely, resulting in exponential decay for 
V and the total center-of-mass kinetic energy. Asymp- 
totically for the relative co-ordinate Eq. (El) becomes 
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This is a system of two coupled non-linear second or- 
der differential equations. Generally such systems do not 
admit closed-form, analytical solution. Somewhat sur- 
prisingly, we now point out that Eq. (^ admits a general 
solution in terms of elliptic functions. 

The approach is standard. First we find two integrals 
of the motion, reducing the four (phase space) dimen- 
sional system in Eq. (|g) to a two dimensional (phase 
space) system. The integrals are the energy and a gen- 
eralization of angular momentum. The inter-ion energy 

results from taking the dot product of Eq. (||) with R, 
forming the total differential, and integrating to find the 
integration constant. 



(6) 



Since the equations have manifest rotational symmetry, 
there is a conserved angular momentum. As always with 
a magnetic field, the total angular momentum receives a 
contribution from the magnetic field. Proceed by taking 
the vector cross product of R and Eq. (0) to find 



dL _ ujQ _ di?^ 
dt 2 ^ ■'' dt 







where, as always, R = \R\, and / = {^)^ is the term due 
to special relativity. The angular momentum per unit 

mass L = z ■ {R X R) ^ R'^^ is the standard definition. 

Now, using the inter-ion energy integral Eq. (^) , / can 
be written entirely as a function of R. Doing so for / in 



(7) 
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Eq. (|^) and integrating leads to the integration constant 



Lo = L 



Wo 



(1 
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Lq represents the generahzed angular momentum. 

Since they are independent, the constants of motion 
in equations Eq. (||) and Eq. (^) constrain the motion 
to lie in a two-dimensional surface in the original four- 
dimensional phase space. Of course, that fact by itself is 
insufficient to guarantee integrability of the equations of 
motion in closed form. However additional peculiarities 
of this system Eq. (||) result in closed form solution. 

In polar co-ordinates the kinetic energy in the potential 
energy equation can be written 
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and solving Eq. (g) for L and substituting we find that 
Eq. (H) becomes. 



dFl -2"° 



4e2 {Lo + aR + PR^ 



where 



mpR 
and /? 



R^ 



(10) 
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Since the 



RHS involves only five consecutive powers of R (namely, 
R^, R, i?", ...R^^). the equation is that of an elliptic func- 
tion. 

More explicitly, we now compute the orbital period of 
the dark state and find the orbit trajectory parametri- 
caly. To compute the period we rewrite Eq. ( p^ ) as 



dt 



dR 



y/u - Ll/R^ - n/R - 2a(3R - /S^R^ 



(11) 



with u — 2uq — 2Lp[3 — o? , and n — 2Lpa + 

The integral is a combination of standard elliptic func- 
tions. In lab co-ordinates i?, the orbits will in general 
be open (with some precession rate which can be written 
in terms of complete elliptic integrals) just as viewing the 
orbits in the i?, t co-ordinates, where now "precession" in 
t in simply the period of the orbit. The period T of these 
orbits is thus given by a contour integral of the RHS of 
Eq. ( pi] ) around the cut running between the classical 
turning points (we label) ap and ai, namely, 



T = \ dt^ 



dR 



yju - LllR^ - n/R - 2af3R - /3^R'^ 



RdR 



iPf ^{R - ao){R - ai){R - a2){R - as) 
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where the are the roots of the fourth degree polyno- 
mial written in Eq. ( p7[ ) . By looking at the signs of terms 
in the polynomial we can see that there can be at most 
two real positive roots. Physically we expect there to be 
exactly two real positive roots which we have called ap 
and ai . These are the classical turning points of the mo- 
tion, and represent the furthest and nearest approaches 
of the particles. 



Furthermore, in the system we are working with, for 
typical values of parameters, we find that all roots are 
real, with two positive and two negative. We may then 
order the roots oq > oi > 02 > 03. Note also that the 
canonical choice of phase for the square root on the cut 
between ap and ai is i and so the period in Eq. ( [l^ ) is 
real and positive. 

Finally, computing the integral in Eq. ( |l^ yields (no- 
tation is that in Ref. jiot). 



2 

r = — 



{ap ~ a3)n(^i ^, A:) -|- aiK{k) 
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where K and H are respectively the complete ellip- 
tic integrals of the first and third kind, and p = 
yj{ap — 02) (ai — aa) and k ~ \/("Q ai^)(a2 03) ^j^^ 
square root of the cross-ratio of the roots. Note that the 
first argument in the H is negative, as it should be on 
physical grounds, since H is convergent for any negative 
argument. 

One of the most striking experimental surprises of the 
two identical ion system is the discovery of cyclotron 
mode- locking^. In these events the two frequency traces 
corresponding (approximately) to the individual ions mo- 
tions meld into one trace. This visible trace is the cen- 
ter of mass motion of the dark state. Our analysis in- 
dicates that there is another invisible (as a dipolc) fre- 
quency branch associated with the inter-ion motion and 
that it has frequency ^ with T of Eq. (|l3|) . For the case 
of two protons in a typical precision Penning trap (at 
LOp ^ 5x10®) we find that Eq. ( |l3|) yields frequencies are 
some tens of Hertz different than up. It would be an inter- 
esting test to apply a sequence of dipole and quadrupolar 
fields to make transitions between dark states and (visi- 
ble) center of mass states. 

By standard means we now derive explicit formulae 
for the shape of the dark state orbits. Recall that, by 
definition of the angular momentum, L, and Eq. (ph 



'^'t> ^ Lp a 
dt^ R-^ R ' 



(14) 



Thus, eliminating time between this and Eq. (jl^) we find 



1 

i/3 



^{R - ap){R - ai)(i? - a2)(i? - as) 



(15) 



which may be evaluated in terms of incomplete elliptic 
functions. We find 



' Lp 
'a2 
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where, again, the arc the (ordered) roots of the poly- 
nomial 

P{R) = -/J^i?^ - 2a/3i?3 + {2up - - 2Lpl3)R^ 
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with a and f3 as defined previously and where 



sin 6l(i?) 
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Note directly from Eq. ( |16| ) and Eq. ( |lg| ) that the pre- 
cession of these orbits is given by twice the RHS Eq. (|l6|) 
with each incomplete elliptic functions replaced by its 
complete elliptic counterpart. 

We have completed a numerical simulation of the sys- 
tem Eq. (0) and Eq. (|) for a range of initial conditions. 
To abet numerical stability those equations were rewrit- 
ten in the co-rotating frame and integrated using com- 
mercial (IDL*'") routines on a DEC Alpha workstation. 
Some of these IDL*™ programs link compiled versions of 
CERN's Mathlib elliptic function routines. The results 
from a typical run are shown in Figures 1 (resp. 2) where 
both the Mo of Eq. (^) (resp. Lq of Eq. (||)) are plotted 
as functions of time. 

The figures show that initially the motions of the ions 
are essentially independent as the energy dissipates. Dur- 
ing this regime the total energy of the system is split 
between the center of mass motion and the inter-ion mo- 
tion. Note that due to the large dynamic range of these 
simulations we have plotted the logarithm of the energy. 
Thus, the linear decay of the envelope of the inter-ion en- 
ergy uq in this initial regime is the exponential damping 
of the energy of the system as a whole. 

Eventually the center of charge motion damps away 
appreciably and the remaining inter- ion motion persists. 
As described earlier, in real experiments of this type the 
dark state we are describing is likely to be effectively the 
final state since we expect the inter-ion motion to decay 
via quadrupole radiation on a timescale long compared 
with typical two-ion experiments. For our simulation this 
final state is reached at simulated time 150, after which 
both Uq and Lq are essentially constant (up to numerical 
accuracy of the simulations). 

In conclusion, we have derived closed form analytic for- 
mulae for the dark state of two identical ions in a Pen- 
ning trap. To find this solution, we assumed that the 
pair is near the center of the trap (we have completely 
neglected the effect of the trap's electrostatic fields) and 
that the motion of the ions is confined to the same az- 
imuthal plane. It is straightforward to include in this 
analysis the effects of the trap's electric field and also a 
fixed average vertical offset between the cyclotron planes 
of the ions. This results in formulae for the two integrals 
of motion that have additional terms compared with the 
Eq. (H) and Eq. (^. However, the resulting equations of 
motion for the dark state are no longer solvable in terms 
of known functions. 
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FIG. 2. The internal angular momentum Lo as a function of time. 



